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Abstract 

The discrete Frechet distance is a useful similarity measure for comparing two sequences of points 
P = (pi,... ,Pm) and Q = (gi,..., g„). In many applications, the quality of the matching can be 
improved if we let Q undergo some transformation relative to P. In this paper we consider the problem 
of finding a translation of Q that brings the discrete Frechet distance between P and Q to a minimum. 
We devise an algorithm that computes the minimum discrete Frechet distance under translation in 
and runs in 0{m^n^{l + log(n/TO)) log(m + n)) time, assuming m < n. This improves a previous 
algorithm of Jiang et al. Col, which runs in 0{m^n^ log(TO + n)) time. 
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1 Introduction 


Shape matching is an important area of computational geometry, that has applications in computer vision, 
pattern recognition, and other fields that are concerned with matching objects by shape similarity. Generally, 
in shape matching we are given two geometric objects A and B and we want to measure to what extent they 
are similar. Usually we may allow certain transformations, like translations, rotations and/or scalings, of 
one object relative to the other, in order to improve the quality of the match. 

In many applications, the input data consists of finite sets of points sampled from the actual objects. To 
measure similarity between the sampled point sets, various distance functions have been used. One popular 
function is the Hausdorff distance that equals to the maximum distance from a point in one set to its nearest 
point in the other set. However, when the objects which we compare are curves, sequences, or contours 
of larger objects, and the sampled points are ordered along the compared contours, the discrete Frechet 
distance may be a more appropriate similarity measure. This is because the discrete Frechet distance takes 
into account the ordering of the points along the contours which the Hausdorff distance ignores. Comparing 
curves and sequences is a major task that arises in computer vision, image processing and bioinformatics 
(e.g., in matching backbone sequences of proteins). 

The discrete Frechet distance between a sequence of points P and another sequence of points Q is de¬ 
fined as fhe minimum, over all possible independenf (forward) fraversals of fhe sequences, of fhe maximum 
disfance befween fhe currenf poinf of P and fhe currenf poinf of Q during fhe fraversals. See below and in 
Secfionl^for a more formal definilion. 

In Ibis work, we focus on fhe problem of compufing fhe minimum discrete Frechef disfance under 
translation. Thaf is, given fwo sequences P and Qofm and n poinfs, respectively, in fhe plane, we wish fo 
franslafe Q by a vecfor t G such fhaf fhe discrefe Frechef disfance befween P and Q + f is minimized. 


Background. The Frechef disfance has been extensively sfudied during fhe pasf 20 years. The main varianf, 
fhe confinuous Frechef disfance, where no fransformafion is allowed, measures similarify befween (polyg¬ 
onal) curves. If is fhe smallesf 6 for which fhere exisf forward simulfaneous fraversals of fhe fwo curves, 
from sfarf fo end, so fhaf af all limes fhe disfance befween fhe corresponding poinfs on fhe curves is al 
mosl 6. The discrefe Frechef disfance considers sequences P and Q of poinfs instead of curves. If is de¬ 
fined analogously, where (a) fhe simulfaneous fraversals of fhe sequences are represenfed as a sequence of 
pairs ■ ■ ■, where G P, G Q, for f = 1,..., /, (b) fhe firsl (resp., lasf) pair 

consisfs of fhe sfarfing (resp., terminal) poinfs of fhe fwo sequences, and (c) each (p(®),q(*)) is oblained 
from q*-®"^^) by moving one (or bolh) poinl(s) fo fhe nexl position in fhe corresponding sequence. 

Mosl sludies of fhe problem consider fhe siluafion where no Iranslafion (or any olher fransformafion) is al¬ 
lowed. In Ihis “sfalionary” case, fhe discrefe Frechef disfance in fhe plane can be computed, using dynamic 
programming, in 0{mn) lime (Filer and Mannila 0)- Agarwal ef al. fH slighfly improve Ihis bound, 

' mn log log n ' 

' time on 


and show fhaf fhe (sfalionary) discrefe Frechef disfance can be compuled in O 


logn 


a word RAM, and a very recenl resull of Bringmann ||5l indicafes fhaf a subsfanlially subquadratic solu¬ 
tion (one fhaf runs in time 0{{mnY~^), for some <5 > 0) is unlikely fo exisf. All and Godau [Jj showed 
fhaf fhe (sfalionary) confinuous Frechef disfance of fwo planar polygonal curves wilh m and n edges, re¬ 
spectively, can be computed, using dynamic programming, in 0{mn\ogmn) time. This has been slightly 
improved recently by Buchin et al. |6|, who showed that the continuous Frechet distance can be computed 
in 0(A^^(log A^)^/^(loglog time on a pointer machine, and in 0(A^^(loglog A^)^) time on a word 

RAM (here N = m = n denotes the number of edges in each curve). In short, the best known algorithms 
for the stationary case, for both discrete and continuous variants, hover around the quadratic time bound. 
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Not surprisingly, the problems become much harder, and their solutions much less efficient, when trans¬ 
lations (or other transformations) are allowed. For the problem of computing the minimum continuous 
Frechet distance under translation, Alt et al. 14] give an algorithm with 0{m^n^{m + n)^ log(m + n)) run¬ 
ning time, where m and n are the number of edges in the curves. They also give a (1 -|- e)-approximation 
algorithm for the problem, that runs in 0{e~‘^mn) time. That is, they compute a translation of one of the 
curves relative to the other, such that the Frechet distance between the resulting curves is at most (l-fe) times 
the minimum Frechet distance under any translation. In three dimensions, Wenk ifldll showed that, given two 
polygonal chains with m and n edges respectively, the minimum continuous Frechet distance between them, 
under any reasonable family of transformations, can be computed in log(m-l-n)) time, where 

/ is the number of degrees of freedom for moving one chain relative to the other. So with translations alone 
(/ = 3), the minimum continuous Frechet distance in can be computed in 0{{m + log(m -|- n)) 
time, and when both translations and rotations are allowed (/ = 6), the corresponding minimum continuous 
Frechet distance can be computed in 0{{m + log(m -|- n)) time. 

The situation with the discrete Frechet distance under translation is somewhat better, albeit still ineffi¬ 
cient. Jiang et al. lITOll show that, given two sequences of points in the plane, the minimum discrete Frechet 
distance between them under translation can be computed in 0{m?n^ log(m + n)) time. For the case where 
both rotations and translations are allowed, they give an algorithm that runs in 0(m^n^ log(m -f n)) time. 
They also design a heuristic method for aligning two sequences of points under translation and rotation in 
three dimensions. Mosig et al. fV2\ present an approximation algorithm that computes the discrete Frechet 
distance under translation, rotation and scaling in the plane, up to a factor close to 2, and runs in 0{m?n?) 
time. 

Our results. Our algorithm improves the bound of Jiang et al. ifT^ by a nearly linear factor, with running 
time + \og{n/m)) log(m -|- n)), assuming m < n. It uses a 0/1-matrix M{P, Q) of size mx n, 

whose rows (resp., columns) correspond to the points of P (resp., of Q). Assuming a stationary situation, 
or, rather, a fixed translation of Q, an entry in the matrix is equal to 1 if and only if the distance between 
the two corresponding points is at most <5, where 6 is some fixed distance threshold. We use (f, j) to denote 
an entry in the matrix that corresponds to the points pi and qj, and we use Mij to denote its value. The 
discrete Frechet distance is at most 6 if and only if there is a row- and column-monotone path of ones in M 
that starts at (1,1) and ends at (m, re) (see Section|^for a more precise definition). 

We can partition the plane of translations into a subdivision As with 0{m?'n?) regions, so that, for all 
translations in the same region, the matrix M is fixed (for the fixed 4). We then traverse the regions of As, 
moving at each step from one region to a neighboring one. Assuming general position, in each step of our 
traversal exactly one entry of M changes from 1 to 0 or vice versa. We present a dynamic data structure 
r(M) that supports an update of an entry of M, in 0(m(l -|- log(re/rre))) time, assuming rre < re|^and then 
re-determines whether there is a monotone path of ones from (1,1) to (rre, re), in 0(1) additional time. If 
we find such a monotone path in M, we have found a translation t (actually a whole region of translation^ 
such that the discrete Frechet distance between P and Q -|- f is at most 5. Otherwise, when we traverse 
the entire As and fail after each update, we conclude that no such translation exists. Using this procedure, 
combined with the parametric searching technique IfTTl . we obtain an algorithm for computing the minimum 
discrete Frechet distance under translation. 

We reduce the dynamic maintenance of M to dynamic maintenance of reachability in a planar graph, 
as edges are inserted and deleted to/from the graph. Specifically, we can think of (the 1-entries of) M as a 
representation of a planar directed graph with N < mn nodes. Each 1-entry of M corresponds to a node 

'This is without loss of generality as we can change the roles of m and n by flipping M. 

^For a critical value of 5, the region can degenerate to a single vertex of As; see Sectionsj^andj^for details. 
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in the graph, and each possible forward move in a joint traversal is represented by an edge (see Section 
for details). Then, determining whether there is a row- and column-monotone path of ones from (1,1) to 
(m, n) corresponds to a reachability query in the graph (from (1,1) to (m, n)). 

A data structure for dynamic maintenance of reachability in directed planar graphs was given by Sub- 
ramanian mi. This data structure supports updates and reachability queries in 0(iV^/^ log A^) time, where 
N is the number of nodes in the graph. Diks and Sankowski |l8l improved this data structure, and gave a 
structure that supports updates and reachability queries in log^ N) time. 

We give a simpler and more efficient structure for maintaining reachability in M that exploits its special 
structure. Our structure can update reachability information in M in 0(m(l -|- log(n/m))) time, assuming 
m < n, and answers reachability query (from (1,1) to [m, re)) in 0(1) time. In contrast, the data structure 
of lUl applied in our context performs an update and a query in 0((rrere )l/2 log^ (rrere)) time. Using our 
structure, we obtain an algorithm for computing the minimum discrete Frechet distance under translation 
that runs in 0(m^re^(l -|- log(re/rre)) log(?re -|- re)) time (again, assuming rre < re). 

To summarize the contributions of this paper are twofold: (a) The reduction of the problem of computing 
the minimum discrete Frechet distance to a dynamic planar directed graph reachability problem, (b) An 
efficient data structure for this reachability problem. For rre « re our structure is faster than the general 
reachability structure of [8] by a polylogarithmic factor, and when rre <C re the improvement is considerably 
more significant (roughly by a factor ^njm). Moreover, our data structure is simpler than that of Diks and 
Sankowski. 

2 Preliminaries 

We now define fhe (stationary) discrete Frechet distance formally. Let P = (pi,...,pm) and Q = 
(gi,..., qri) be the two planar sequences defined in the introduction. 

For some fixed disfance J > 0 we define a 0/1-mafrix Ms{P, Q) formally as follows. The rows (resp., 
columns) of Ms{P, Q) correspond to the points of P (resp., of Q) in their given order. An entry {i,j) of 
Ms{P, Q) is 1 if the distance between pi and qj is at most 6, and is 0 otherwise, we denote Ms{P, Q) by M 
when P and Q and 6 are clear from the context. 

The directed graph SsiP, Q) associated with P, Q and <5 has a vertex for each pair {pi,qj) G P x Q and 
an edge for each pair of adjacent ones in M^{P, Q). Specifically, we have an edge from {pi, qj) to {pi+i, qj) 
if and only if both (z,y) and (i -|- 1, j) are 1 in M, an edge from {pi, qj) to {pi, qj+i) if and only if both {i,j) 
and (i, j -|- 1) are 1 in M, and an edge from {pi, qj) to (pi+i, qj+i) if and only if both {i,j) and {i + lj + l) 
are 1 in M. we denote Gs{P, Q) by G when P and Q and 5 are clear from the context. 

The (stationary) discrete Frechet distance between P and Q, denoted by (5* (P, Q), is the smallest <5 > 0 
for which {pm,qn) is reachable from {pi,qi) in Informally, think of P and Q as two sequences of 
stepping stones and of two frogs, the P-frog and the Q-frog, where the P-frog has to visit all the P-stones 
in order and the Q-frog has to visit all the Q-stones in order. The frogs are connected by a rope of length 
6, and are initially placed at pi and qi, respectively. At each move, either one of the frogs jumps from its 
current stone to the next one and the other stays at its current stone, or both of them jump simultaneously 
from their current stones to the next ones. Furthermore, such a jump is allowed only if the distances between 
the two frogs before and after the jump are both at most 6. Then (5* (P, Q) is the smallest <5 > 0 for which 
there exists a sequence of jumps that gets the frogs to pm and q^, respectively. 

The problem of computing the minimum discrete Frechet distance under translation, as reviewed in the 
introduction, is to find a translation t such that 5*(P, Q + f) is minimized. 

We say that an entry (i,j) of M is reachable from an entry {k,l), with k < i,l < j, if {pi,qj) is 
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reachable from in S- A path from {p^, qi) to {pi, qj) in S corresponds to a (weakly) row-monotone 

and column-monotone sequence of ones in M connecting the one in entry (fc, 1) to the one in entry (i, j). 
This is sequence consists of three kinds of moves: 1) upward moves between entries of the form (r, s) to 
(r -|- 1, s) in which the P-frog moves from pr to pr+i, 2) right moves between entries of the form (r, s) to 
(r, s -|- 1) in which the P-frog moves from qg to qs-i-i> and 3) diagonal moves between entries of the form 
(r, s) to (r -|- 1, s -|- 1) in which the P-frog moves from qg to both frogs move simultaneously — the 
P-frog from pr to Pr+i, and the Q-frog from qg to qs+i- See Figure[T] We call such a monotone sequence 
of ones in M a path in M from {k, 1) to (i, j). To determine whether d*{P, Q) < 6, we need to determine 
whether there is such a path in M that starts at (1,1) and ends at {m,n). We say that an entry (i, j) of M is 
reachable if there is a path from (1,1) to (i, j). 

We denote the concatenation of two paths tti, 7 r2 by tti • 7 r2, assuming that the last entry of vri is the first 
entry of 7r2; this entry appears only once in the concatenation. 
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Figure 1: (a) A reachability path from (4, 3) to (7, 7). (b) A reachability path from (1,1) to (8,12). 


We use a decomposition of M into blocks. A block is a submatrix of M that corresponds to contiguous 
subsequences P' and Q' of P and Q, respectively. We denote by M{P', Q') the block of M(P, Q) formed 
by P' and Q'. Consider a block M{P', Q') of M(P, Q). Let p~ (resp., p'^) denote the first (resp., last) 
point of P' and let q~ (resp., 9 +) denote the first (resp., last) point of Q'. We call the entries of M{P', Q') 
corresponding to {p~} x Q' U P' x {<?“} the input boundary of M{P', Q'), and denote it by M{P', Q')~ 
(the common entry corresponding to {p~} x {q“} appears only once in the boundary). Similarly, we call the 
entries of M (P', Q') corresponding to {p"''} x Q'uP' x {q"''} the output boundary of M (P', Q'), and denote 
it by M{P', Q')^ (with a similar suppression of the duplication of the common element {p+ } X {q+}). 
Note that there is a two-entry overlap between the input and output boundaries. We enumerate the entries 
of M(P', Q')~ by first enumerating the entries of {p~} x Q' from right to left (i.e., backwards) and then 
the remaining entries of P' x {q~ } from bottom to top (forward). We enumerate the entries of M (P', Q')'^ 
by first enumerating the entries of P' x { 9 “'"} from bottom to top (forward) and then the remaining entries 
of {p}+ X Q' from right to left (backwards). Informally, M{P', Q') is enumerated in “clockwise” order, 
while M{P' is enumerated in “counterclockwise” order; see Figure]^ For two entries i,j of this 
enumeration of an input or output boundary B of M (P', Q'), we use [i, j] to denote the sequence of entries 
(i,i-hl,...,j)ofP. 

We also use the following definitions. We call the entries corresponding to P' x {g“} the vertical 
input boundary of M{P', Q'), and denote it by M{P', Q') ~. We call the entries corresponding to P' x 
{(?■'■} the vertical output boundary of M{P', Q'), and denote it by M{P', Q') +. That is, M{P', Q') ~ and 
M{P', Q') + are the vertical parts of M{P\ Q')~ and M{P', respectively. We enumerate the entries 
of each vertical boundary from bottom to top. 
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Figure 2: (a) A (highlighted) block M(P', Q') of M(P, Q), where P' = {p 5 ,P 6 ,P 7 ,P 8 ) and Q' = (gy, 98 ,, 912). 
(b) The input boundary M{P', Q')~ and the output boundary M{P', Q')'^ of M{P', Q') are marked, with the order¬ 
ings of their elements. 


3 The subdivision of the plane of translations 

We first consider the corresponding decision problem. That is, given a value <5 > 0, we wish to decide 
whether there exists a translation t such that 6*{P, Q + t) < 6. 

For a point X G let Ds{x) be the disk of radius S centered at x. Given two points pi G P and 
Qj G Q, consider the disk Ds{pi — Qj), and notice that t G Ds{pi — qj) if and only if \\{pi — qj) — t\\ < <5 (or 
\\Pi ~ {Qj + ^) II < <5)- That is, Dsipi — qj) is precisely the set of translations t for which qj + f is at distance 
at most 6 from pi. 

We construct the arrangement As = As{P, Q) of the disks mV = {Ds{pi — qj) \ {pi, qj) G P x Q}. 
We assume general position of the points. That is, we assume that (a) no more than two boundaries of these 
disks intersect in a common vertex of As, and (b) no pair of the disks are tangent to each other. Nevertheless, 
such a degeneracy can arise when 5 is a critical value (see Section for details about critical values of 5 
that arise during the optimization procedure), but we assume that at most one such degeneracy can happen 
for a given 5. Since the number of disks is mn, the combinatorial complexity of As is Let / be 

a face of As of any dimension 0,1 or 2 (by convention, / is assumed to be relatively open), and let t e / be 
a translation. Then, for points pi G P, qj G Q, qj + f is at distance at most S from pi if and only if the disk 
Ds{pi — qj) contains / (otherwise, the disk is disjoint from /). Since this holds for every f G /, it follows 
that / corresponds to a unique pairwise-distances matrix M{P, Q + t), for any t G f. We denote this matrix 
by M{P, Q + /), for short. 

The setup just described leads to the following naive solution for the decision problem. Construct 
the arrangement As for the given distance <5, and traverse its faces. For each face f G As, form the 
corresponding pairwise-distances matrix M{P, Q + /), and solve the (stationary) discrete Frechet distance 
decision problem for P and Q + f using a straightforward dynamic programming on M{P, Q + f) (or the 
more sophisticated slightly subquadratic algorithm of Agarwal et al. ||T1). If <5* {P, Q + f) < 6 for some face 
/, we conclude that there exists a translation t such that <5* {P, Q -\-t) < 6 (any translation t G f would do). 
If the entire arrangement As is traversed and no face / of As satisfies 5* (P, Q + f) < 6,^/0 determine that 
6 *{P, Q + t) > (5 for all translations t G M^. The complexity of As is 0{m?n^), and solving the discrete 
Frechet distance decision problem for each face of As takes 0{mn) time (or slightly less, as in |Tj). Hence, 
the solution just described for the decision problem takes (slightly less than) O(m^n^) time. 

Jiang et al. ifTOll used an equivalent solution for the decision problem, that takes the same asymptotic 
running time. Rephrasing their procedure in terms of As, they test whether (i* (P, Q +1) < 5 for translations 
t corresponding to the vertices of As, and over an additional set of mn translations, one chosen from the 
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boundary of each disk. The correctness of this approach follows by observing that if / is a face of As and t 
is any point on df, then all the 1-entries of M (P, Q + f) are also 1-entries of M (P, Q -h t), so it suffices to 
test the vertices of /, or, if / has no vertices, to test an arbitrary point t G df. We will use this observation 
in our implementation of the optimization procedure. 

Our naive solution is similar to the algorithm of Jiang et al. ifTOl . in the sense that they both discretize 
the set of possible translations. However, our solution is more suitable for the improvement of this naive 
bound, that we present in Section]^ since it allows us to traverse the set of possible translations in a manner 
that introduces only a single change in M{P, Q + /), when we move from one face / of translations to a 
neighboring one. 

To exploit this property we need a data structure that maintains reachability data for M, and updates 
it efficiently after each change. We present this structure in two stages. First, in Section we present a 
compact reachability structure for blocks of M(P, Q + /), which is the main building block of the overall 
structure. Then, in Section we present the overall data structure, and show how to use it to improve the 
naive solution sketched above by a nearly linear factor. 


4 Compact representation of reachability in a block 


Let P be a block of M = M(P, Q + f) of size r x c, and suppose that we have already computed the 
reachable entries of B~ and we then wish to compute the reachable entries of P+. If the entries of the block 
are given explicitly, this can be done in 0{rc) time using dynamic programming (or slightly faster using the 
algorithm of |[T]). Our goal in this section is to design a data structure, that we denote as <F(P), that allows 
us to compute the reachable entries of from the reachable entries of B~, in 0(r -|- c) time. The overall 
data structure itself is constructed recursively from these block structures (see Section]^ for details), and 
implicitly accesses all the entries of B. The advantage of using this block decomposition is that updating 
the structure can be done more efficiently. 
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Figure 3: (a) Two entries i, j of B~ and two entries a{i) and a{j) of that are reachable from i and j, respectively. 
Since i < j and a{i) > cr(j), cr{j) is reachable also from i, and a{i) is reachable also from j. (b) The intervals 
[aA{k),(Jz{k)], for any 1-entry k of B~, are either disjoint or overlap in a common subinterval. Neither of these 
intervals can sffictly contain both endpoints of the other. These intervals are defined (and shown in the figure) only for 
1-entries of B~. 
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Observation 4.1. Let B be a block ofM and let i,j be two entries ofB~ such that j > i (in the “clockwise ” 
order defined in Section^. Let a{i) be an entry of B^ that is reachable from i, and let a(j) be an entry of 
B^ that is reachable from j. If cr[j) < cj(i) (in the corresponding “counterclockwise” order) then a(j) is 
also reachable from i, and a{i) is also reachable from j. 

Proof See Figurej^a). Since a{i) is reachable from i, there is a (monotone) path 7r(i, erf)) from i to a{i) in 
M. Similarly, since a{j) is reachable fromy, there is a (monotone) path 7r(j, cr{j)) fromy to a{j). Since i < 
j and a{j) < a(i), 7r(i, cr(i)) must cross Tr(j, cr(j)) (i.e., there exists a 1-entry e G irf, a{i)) n 7r(j, cr{j)). 
Hence, 7r(z, cr{i)) can be decomposed into two subpaths 7r(i, e), 7r(e, a{i)) such that 7r(z, af)) = 7r(z, e) • 
7r(e, erf)), and 7r(y, cr{j)) can be similarly decomposed as 7r(j, cr{j)) = 7r(j, e) • 7r(e, er(j)). As a result, the 
paths 7r(i, cr{j)) = 7r(i, e) • 7r(e, cr{j)) and 7r{j, a(i)) = 7r(j, e) ■ 7r(e, cr(i)) are also (monotone) paths, and 
the claim follows. □ 

Corollary 4.2. Let B be a block of M, let i be an entry of B~ and let eri(i)^a 2 {i) be two entries in B^ 
that are both reachable from i, with cri{i) < a 2 {i)- If there exists an entry a{j) that is reachable from some 
j G B~, such that eri{i) < cr(j) < cr 2 (i), then cr(j) is also reachable from i. 

Proof By Observation |4.1[ if i < j, then a{j) is reachable from i since a{j) < cr 2 {i). If i > j, then a{j) 
is reachable from i since cji (i) < a{j). □ 

The corollary is applied as follows. Let i be an entry of B~, let cr^(i) and azf) denote the first and 
last entries in B~^ that are reachable from i. (Note that for these entries to be defined, fhe value of the entry 
i must be 1. Symmetrically, the values of both aAf) and azf), if defined, must be equal to 1.) Then the 
interval [aAf), erzf)] can only contain entries of the following three types. 

1. 1 -entries that are reachable from i. 

2. 0-entries. 

3. 1-entries that are not reachable from i, nor from any other entry of i?“. 

In other words, [aAf), erz{i)] cannot contain 1-entries that are reachable from some j in B~ and not 
from i. 

Corollary 4.3. Let B be a block of M and let i and j be two entries of B~ such that j > i. Then 
ctaU) > CFAf) and az{j) > crz{i)- 

Proof Assume to the contrary that <T^(f) > cr^(i). Then, according to Observation |4T| (Ta{3) is reachable 
from i. Hence, crAf) < ctaU), a contradiction. Similarly, if az{j) < erzf), then Observation |4. 1 1 implies 
that azf) is reachable from j. Hence, <7z{j) > erzf), contradiction. □ 

In other words, [aAi^j^zf)] and [crAij),o'z(j)] can be either disjoint or overlap in a common subinter¬ 
val, but they cannot be properly nested inside one another (that is, neither of these intervals can contain both 
endpoints of the other in its “interior”). Note, however, that one interval can weakly contain the other. That 
is, if one interval contains the other, then either aA{i) = ^aU) or az{i) = crz{j), or both. See Figurej^b). 

Let H be a block of M of size r x c. We construct a data structure 4>(H) for B, which stores the 
following information. (Here we only specify the structure; its construction is detailed in Section]^) 

1. For each 1-entry i of B~ we store 

(a) the first entry aAf) of B~^ that is reachable from i, and 
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(b) the last entry azii) of reachable from i. 

2. For each 1-entry j of B~^ we store 

(a) a flag /(j) indicating whether j is reachable from some entry of B~. 

(b) a list LaU) of the 1-entries i G B~ such that aA{i) = j, and 

(c) a list Lz{j) of the 1 -entries i G B~ such that az{i) = j- 

Lemma 4.4. Given the data structure ^{B)for a block B, and given the entries of B~ that are reachable 
from ( 1 , 1 ), we can determine, in 0{r + c) time, the entries of B^ that are reachable from ( 1 , 1 ). 


Proof We go over the reachable 1-entries of B~ in order. For each such entry i, we go over the entries in the 
interval I{i) = [max{crz(z“), (T^(i)}, (Tz{i)] of B^, where i~ is the previous reachable 1-entry of B~ (for 
the first reachable entry i of B~, I{i) = [aA{i)-,(^z{i)] and i~ is undefined). Note that, by Corollary 4.3 


max{az{i~),crA{i)} G [aA{i),crzii)] so I(i) C [aA(i),o'zii)]- (The entries of [cr^(i), (Tz(i)] that precede 
m.ax.{(Tz{i~),(TA{i)'\ were already processed when we went over I{i~) or over intervals associated with 
earlier indices.) 

For each 1-entry j of I{i) that is reachable from some entry of B~ (according to the flag f{j)), we 
determine that j is reachable also from (1,1). Since we traverse each interval [cTA(i), c’'z(*)] starting from 
max{crz(i~), (Tyi(i)}, the internal portions of the subintervals that we inspect are pairwise disjoint, implying 
that the running time is linear in r -|- c. We omit the straightforward proof of correctness of this procedure. 

□ 


5 Dynamic maintenance of reachability in M{P,Q + /) 


We present a data structure, that uses the compact representation of reachability in a block of the previous 
section, to support an update of a single entry of M in 0(m(l -|- log(n/m))) time, assuming m < n. We 
present this data structure in two stages. First, in Section [5TT] we show how to support an update of a single 
entry in 0{m) time, in the case where M is a square matrix of size m x m. Then, in Section 5.2 we 


generalize this data structure to support an update of a single entry in 0 (m(l -|- log(n/m))) time, in the 
general case where M is an m x n matrix with m < n (the case m > n is treated in a fully symmetric 
manner). 

In Section |5.3[ we describe the overall decision procedure that improves the naive solution sketched in 
Section]^ using this dynamic data structure. 


5.1 A dynamic data structure for reachability maintenance in a square matrix 

We store the reachability data of M (P, Q + f) (of some arbitrary face / from which we start the traversal 
of the arrangement Af) in a so-called decomposition tree F, by halving P and Q alternately. That is, the 
root u of F corresponds to the entire matrix M{P,Q + f) and we store at v the reachability information 
<h(M(P, Q -|- /)), as described in the previous section. (The actual construction of the reachability data, 
at all nodes of F, is done bottom-up, as described below.) In the next level of F we partition P into two 
subsequences Pi,P2, of at most [m,/ 2 j -|- 1 points each, such that the last point of Pi is the first point 
of P2, and obtain a corresponding “horizontal” partition of M(P, Q + f) into two blocks M(Pi, Q + /), 
M{P 2 , Q + f), each of size at most ([m/ 2 j -|- 1 ) x m, with a common “horizontal” boundary. We create 
two children vi,V2 of v and store at each Vi the reachability information <I>(M(Pi, Q + /)), for i = 1, 2. In 




By = M{P,UP^,Qy) 



(a) (b) 


Figure 4: A block By = M{Py, Qy), corresponding to a node y of F, is composed of the blocks of the children 
V, w of y. The block By = M{Py, Qy) corresponding to the left child v lies below the block By, = M{Py,, Qy,) that 
corresponds to the right child w, and we have Py = PyU Py, and Qy = Qy = Qy,. 

(a) cr^(z) and a^(k) are examples of reachable entries of B^ and we have a'^(i) = a'^{a'"{i)) and a^(k) = a'"{k). 
cr'“{j) is an entry of that is reachable from B~, but it is not a reachable entry of (from By) since all the paths 
in By that lead to cr'^{j) go through entries of B^ that are not reachable from B~. 

(b) a\{e) = cr^(cr^(f)) and cr|.(£) = cr“(cr|,(£)). 


the next level of F, we partition Q into two subsequences Qi,Q 2 , of at most [m/2\ + 1 points each, such 
that the last point of Qi is the first point of Q 2 , and obtain a corresponding “vertical” partition of each block 
M{Pi, Q + f),i e {1,2}, into two blocks M{Pi, Qj + f),j G {1, 2}, each of size at most ([m/2j +1) x 
([m/2j +1), with a common vertical boundary. We construct four respective grandchildren, and store the 
corresponding reachability structures 4>(M {Pi, Qj + /)) at these nodes. We continue recursively to partition 
each block by halving it horizontally or vertically, alternately, in the same manner, until we reach blocks of 
size 2x2. For each node v of F, let Py and Qy denote the subsequences of P and Q that form the block 
M{Py, Qy + /) that is associated with v. To simplify the notation, we denote 4>(M {Py, Qv + f)) as 4>^, for 
each node v. 

The reachability data ^y at the nodes u of F is computed by a bottom-up traversal of F, starting from the 
leaves. The construction of ^{M{Py,Qy + /)) at a leaf v is trivial, and takes constant time. The following 
lemma provides an efficient procedure for constructing the reachability data at inner nodes of F. 

Lemma 5.1. Let y be an inner node ofT with left and right children v and w, where the blocks stored at v, w 
have a common horizontal boundary. Given the reachability data the data can be computed in 

0{\Py\ + \Qy\) time. An analogous statement holds when the common boundary of the children blocks is 
vertical. 

Proof. Note that in the setup of the lemma, we have Qy = Qy = Qw and Py = Py U Py,. By construction, 
M{Py, Qy) lies below M{Py,, Qy). Denote M(P^, Qy) by By, M{Py,, Qy) by By,, and M{Py, Qy) by By. 
For each entry i of Bf, denote by a\{i) (resp., a^{i)) the first (resp., last) entry of that is reachable 
from i. We also use a^{i) to denote an entry of Bi^ that is reachable from i in By. Analogous notations are 
used for the children blocks By, By,. See Figure]^ 

We first copy the reachability information from the boundaries of By and By, to the boundary of By 
(except for the “interior” portion B*y, of the common boundary Bf n Bf of By and By,, which is not a 
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boundary of By). The data for the 1-entries on the left boundary of Byj (which are of type[T]in the definition 
of <h) is still valid, since the reachability paths of By that start at these entries are fully contained in Bw 
Similarly, the data for the 1-entries on the right boundary of By (which are of type is still valid, since 
the reachability paths of By that end at these entries are fully contained in By. We thus need to determine 
the reachability information from the 1-entries of the input boundary By of By to the entries of the output 
boundary B^ of By,, and merge it with the already available data, to get the complete structure at y. 

First note that an entry j of B^ that is reachable from B~ may now become unreachable from By . 
This happens if all the reachability paths in By to j go through entries on that are not reachable from 
By . See Figure l^a). We thus need to turn the flag /(j) of such enfries fo false. To do fhis, we go over fhe 
enfries of B^ in order, and mainfain a queue Q fhaf safishes fhe invarianf thaf, when we are af an enfry j of 
B^, Q confains all fhe enfries i of B~ fhaf are reachable from By, such fhaf j is reachable from i. Thaf is, 
Q confains all fhe enfries i G Byy, fhaf are reachable from By such fhaf j is reachable from i, and all fhe 
enfries i G B~ \ Byy, (fhaf is, fhe leff side of B~) such thaf j is reachable from i. We sfarf wifh an empfy 
queue. For each l-enfry j of S+ we hrsf go over fhe lisf LaU) (of and for each elemenf i in LaU) 
fhaf is in Byy,, we check if if is reachable from B~ (using fhe flag f{i) from <1>„). If if is, we puf if in Q. 
We also add fo Q each elemenf in LaU) thaf is in B~ \ Byy,. If Q is empfy, fhere is no reachabilify pafh 
from By fo j and we sef f{j) fo be false. We fhen go over fhe lisf Lz{j) (of and remove from Q each 
elemenf in Lz{j) thaf is in Q. This fraversal fakes 0{\Py \ + \Qy\) fime, since each elemenf of B~ appears 
af mosf once in fhe lisfs La and af mosf once in fhe lisfs Lz. The correcfness follows from fhe invarianf fhat 
when we go over an enfry j G B^, all fhe enfries of B~ fhaf j is reachable from, and fhaf are reachable 
from By, are in Q. The invarianf is mainfained correcfly because each fime fhaf an inferval [(TA(f), o'z(f)] 
of an enfry i G B~ begins (and i is reachable from By), i is inserted info Q, and when fhe inferval ends, i 
is removed from Q, so f is in Q for all enfries j fhaf are reachable from i. In conclusion, if Q is empfy, j is 
nof reachable from By and fhe flag f(j) can be fumed false. Ofherwise, i is reachable from By . 

We now update fhe intervals [aA{i),crz{'i)] of the entries i G By and, in correspondence, the lists 
L^((T(f)), L^(cr(z)) of B^ (where a{i) is any entry of B^ that is reachable from i). Consider a 1-entry i 
of By and consider an entry (T'"{i) in Byy,-, that is, cr'^{i) is a 1-entry in [(T^(i), that is reachable from 

i. By transitivity, the entries a‘^{a'"{i)) of that are reachable from cr’'(f) are also reachable from i. We 
update [aA{i)i c’'z(f)] according to this rule, as follows (see Figure^b)). We set (T\{i) = a'^{a\{i)), for 
each entry i G By such that (j\{i) G Byy,-, correspondingly, we also add i to LA(cT^(i)). Similarly, for 
each entry i G By such that a'^{i) G Byy,, we set a^ii) = and we add i to Lz{o-^{i)). (Recall 

that if (T^(f) (or a^{i)) is in B:^ \ Byy,, this reachability information was already copied to ^y and that the 
reachability information for B~ \ Byy, was also copied to <l>y.) Clearly, for each entry i G By, no entry of 
[(T^(i), cj^(i)] is reachable from i. This traversal takes 0(|Pj/| -|- \Qy\) time. 

Finally, when we copied information from B^ to B^, we also copied the lists La and Lz that may 
include entries of Byy,. Since Byy, is not a part of the boundary of By, we need to remove this information 
from the lists La and Lz of B^ . We thus go over the entries of Byy,. For each entry e of B^y,, we remove 
e from L^(fT^(e)) and from Lzicr^ie)). Clearly, this traversal takes 0{\Qy\) time. □ 


We now show how to use Lemma 5.1 to construct F in 0{m?‘) time and to update it, when a single entry 
changes, in 0{m) time. We also show how to determine, using F, whether (m, m) is reachable from (1,1) 
in constant time after the update. 


Lemma 5.2. (a) Given a square matrix M, the decomposition tree F (including the reachability data at its 
nodes) can be constructed from scratch in 0{m‘^) time, (b) If a single entry of M is updated, then F can be 
updated in 0{m) time, (c) Given F, we can determine whether {m, m) is reachable from (1,1) in constant 
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time. 


Proof, (a) We construct F in a bottom-up manner, as prescribed in Lemma 5.1 For the blocks at the leaves, 
the reachability data is computed in brute force, in 0(1) time per block, and at each inner node y, the data is 
computed from the data at its children in time 0{\Py \ + \ Qy\), using Lemma 5.1 we refer to \Py \ + \Qy 


as 


m\ 


the size of the block By at y. The sizes of the blocks at levels 2j — 1 and 2j is O j, and the number of 

these blocks is 0(2^-^). The height of T is [log m]. The cost of the overall construction of T is proportional 
to the sum of the sizes of its blocks (this also holds at the leaf level), which is thus 



: O (m^) . 


(b) The main observation here is that to update T when an entry e of M changes, it suffices to update the 
reachability data along the single path of T of those nodes y for which e £ By. (Actually, because of the 
overlap between block boundaries, there are two such paths that meet at the unique node y for which e 
belongs to the “interior” of the common boundary of the blocks of its children.) The reachability data of the 


nodes along this path is constructed again from scratch in a bottom-up manner, using Lemma 5.1 The cost 
of the updates of these blocks is proportional to the sum of their sizes, which is 


= 0{m). 


(c) To determine whether (m, m) is reachable from (1,1), we simply check in the reachability data structure 
<h(M) of the root of T whether (m,m) is a 1-entry that belongs to [f7A((l, 1 )),<t^((1, 1))] and the flag 
/((m, m)) is true. □ 



5.2 A generalized structure for arbitrary matrices 

We next describe a modified variant of the structure for the case where m and n are unequal. In what follows 
we assume, as above and without loss of generality, that m < n. 

We first partition M into k = Oinjm) square blocks Bi,B 2 ,... ,B}^, of size m x m each such that 
consecutive blocks overlap in a single column. (The last block may be of smaller width, but we handle it 


in the same manner as the other blocks; it is easy to show that the bounds of Lemma 5.2 still hold.) We 


build the decomposition tree and the associated reachability data for each of these blocks, as in Section [5T] 
denote the structure for block by Tj, for z = 1,..., A:. 

We now combine the structures fi,..., f^ into a single global structure f. For this, we construct a 
balanced binary tree T, with k leaves vi,... ,Vk, where Vi, for i = I,... ,k, corresponds to Bi and stores 
fj. Each node v of T represents a block By that is the concatenation of the blocks stored at the leaves of 
the subtree rooted at v. Since each leaf block spans all the rows of M, the common boundary of any pair of 
consecutive blocks consists only of a full single column of M. The same holds at any node y of T, with left 
child V and right-child w. That is, the common boundary Byy, := Bf n B~ between By and By, is vertical, 
and consists of a full single column of M. 

We claim that we can merge the reachability structures of By and ^y, of By, into the structure ^y of 
By in 0{m) time, instead of 0{\By\) = 0{m + \Qy\) time (as was the cost in the preceding subsection), 
which can be much larger. The main observation that facilitates this improvement is that there is no need to 
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maintain the reachability data at the horizontal portions of the boundary of any of the blocks By. This 
follows from the obvious property that any path vr from the initial entry (1,1) to any entry {i,j) in any leaf 
block reaches {i,j) by crossing all the vertical boundaries B 12 , B 23 , ■ ■ ■ that delimit all the preceding leaf 
blocks, and the portion iri of vr within each of the preceding blocks Bi connects an entry on the left vertical 
boundary of Bi to an entry on its right vertical boundary. Note that tti can “crawl” along the lower or upper 
boundary of Bi, but to exit Bi it has to cross the vertical boundary, possibly through its entries in row 1 or 
row n. 

Figure|^is an illustration of an inner block By of f that is composed of a left block By and a right block 

By,. 


By = M{P, Qy U Qy,) 

By = AI{P., Qy) B„ = M{P. Qy,) 



= B* =B- 

Figure 5: A block By = M{Py, Qy), corresponding to a node y of f, is composed of the blocks By = M{Py, Qy) 
and By, = M{Py,, Qy,) of the children v,w of y, with v being the left child and w being the right child. We have 
Py — Py — Pw — B and Qy — Qv U Qy,. 

We therefore use the same reachability data structure at v as defined in the previous subsection, 
except that we limit the input and output domains of its maps to the vertical boundaries only. Recall our 
notation from Section]^ where the left (resp., right) vertical boundary of a block B is denoted as B~ (resp., 
B^). Specifically, denoting the modified structure as = ^(By), it stores the following items. 

1. For each 1-entry i of B~ we store 

(a) the first entry aA{i) of By that is reachable from i, and 

(b) the last entry az{i) of B)^ that is reachable from i. 

2. For each 1-entry j of B)^ we store 

(a) a flag /(j) indicating whether j is reachable from some entry of B~. 

(b) a list LaU) of the 1-entries i G B~ such that d'yi(f) = j, and 

(c) a list Lz{j) of the 1-entries i G B~ such that az{i) = j. 

In other words, ^{By) is a constrained variant of ^{By), obtained by replacing By and by By and B)^, 
respectively. The structure of the root u of a child Fj of f is obtained from <1>„ by first setting, for each 
entry i of By for which cJA(f) is in and az{i) is not in B)^, az{i) to be the last reachable entry k of 
.B+, then updating Lz{k) accordingly, and finally ignoring the horizontal parts of the boundaries of By and 
deleting the data regarding them from the lists La and Lz of entries of B)^. 
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We next claim that the modified structures are sufficient for obtaining reachability data for the blocks 
of f, in the precise sense stated below, and that the structure at an inner node y of f can be obtained from 


the structures at the children of y in 0{m) time. Concretely, we have the following variants of Lemmas 4.4 
and 15. II 


Lemma 5.3. Given the data structure ^{B) for a block B of size r x c., and given the entries of B that 
are reachable from (1,1), we can determine, in 0{r) time, the entries of B^ that are reachable from (1,1). 

Lemma 5.4. Let y be an inner node of f with left and right children v and w. Given the reachability data 
Qv)), Qw)), the data ^{M{Py, Qy)) can be computed in 0{m) time. 

The proof is essentially identical to those of Lemma |4.4| and Lemma |5.1[ except that we restrict the 
domains and the images of each of the maps (i.e., aA,crz, f, La, and Lz) to the vertical portions of the 
boundaries. This is justified using the observation made earlier that all the reachability paths traverse only 
vertical boundaries of the relevant blocks — those that are stored at f, from its leaves up, which span the 
entire range P of rows of M. Since we only traverse vertical boundaries, the cost of constructing ^y from 
and is 0{m). □ 

The following lemma extends Lemma [5^ 

Lemma 5.5. (a) Given the matrix M, T can be constructed in 0{mn) time, (b) If a single entry of M is 
updated, then T can be updated in 0{m{l + log(n/m))) time, assuming m < n. (c) Given T, we can 
determine whether (m, n) is reachable from (1,1) in constant time. 


Proof, (a) We construct the structure fj for each block Bi, and extract from it. We then construct ^y 
for each inner node y £ Tby merging the corresponding data structures of the children of y in 0{m) time. 
We obtain T at the root of T. Since T is of size 0{n/m) and we spend 0{m) time at each block, it takes 
0(n) time to construct f from the leaf structures f* for i = 1,..., /c. The cost of constructing each f j, 
i = 1,... ,k = [n/m], is O(m^), by Lemma 5.2 for a total of 0{krn?‘) = 0{mn). It follows that that the 
overall construction of f takes 0{n + mn) = 0{mn) time. 

(b) To update f when an entry e of M changes we first need to update the reachability structures along a path 
in the structure of the block Bi containing e (If e is in the common column of two blocks we update both 
structures.). This takes 0{m) time by Lemma 5.2 Once we have the updated f j we update the reachability 
structures along the path tt of T of those nodes y for which e £ By. (There are two such paths if e is in 
the common column of two consecutive blocks.) Since the depth of T is 0(log(n/m)) and we spend 0{m) 
time to reconstruct the structure at each node of T, we update T in 0(m(l + log(n/m))) time. 

(c) As in Lemmato determine whether (m, n) is reachable from (1,1), we simply check in the reacha¬ 

bility data structure <I>(M) of the root off whether (m, n) £ [d'yi((l, 1)), dz{{l, 1))] and the flag f{{m, n)) 
is true. □ 


5.3 The overall decision procedure 

We now put together the pieces of the decision procedure. We construct the arrangement .4^ of the disks 
Dsipi — Qj) as in Section]^ in log(m + n)) time. We pick an arbitrary (0-, 1-, or 2- dimensional) 

face /o of As. /o corresponds to a unique matrix M{P, Q + /o) and we construct the data structure f of 
Section 5.2 based on M (P, Q + /o). We then perform a traversal of the entire arrangement As. In each step 
of the traversal we move from a face / of As to a neighbor face f (both faces are of any dimension 0 , 1 , or 
2). In this step, we either enter a single disk of As or exit a single disk of As. This corresponds to a change 
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in a single entry of M(P, Q + /). We update f accordingly, in time 0(m(l + log(n/m))), and thereby 
determine whether (5*(P, <3 + /O < We continue in this manner till we process the entire arrangement. If 
we encounter a face / along the traversal at which (5*(P, Q + /) < we report that the minimum distance 
under translation is < 5, and otherwise we report that the minimum distance is > b. We thus obtain the 
following intermediate result. 

Theorem 5.6. Let P, Q be two sequences of points in of sizes m and n, respectively and let 6 > 0 
be a parameter. Then the decision problem, where we want to determine whether there exists a translation 
t G such that 6*{P, Q + t) <6, can be solved in + log(n/m))) time, assuming that m < n. 

6 The optimization procedure 

We now show how to use the decision procedure of Section to compute the minimum discrete Frechet 
distance under translation. Assume without loss of generality that m < n. As we increase <5, the disks 
Ds {pi — qj) expand, and their arrangement As varies accordingly. Nevertheless, except for a discrete set of 
critical values of b, the combinatorial structure of As does not change. That is, the pairs of intersecting disk 
boundaries remain the same, all their intersection points remain distinct and vary continuously, and no pair 
of disks are tangent to each other. Consequently, the representation of As that we use, namely, a collection of 
circular sequences of vertices, each containing the vertices of As along some circle Cs{pi — qj), for pi G P, 
Qj G Q, sorted along the circle, remain unchanged. The critical values of b, at which this representation of 
. 4.5 changes qualitatively, are 

1. The radii of the disks that have three points of P — Q = {pi — qj | pj G P, qj G Q} on their 
boundaries. 

2. The half-distances between pairs of points of P — Q. 

There are 0{mfiA) critical values (most of which are of type[^, so we cannot afford to enumerate them 
and run an explicit binary search to locate the optimal value of <5 among them. 

Instead, we use the parametric searching technique of lITTIl . In general, using parametric searching can 
be fairly complicated, since it is based on a simulation of a parallel version of the algorithm. However, 
we only have to simulate, by a parallel algorithm, the part of the decision procedure that depends on (the 
unknown value of) the optimum bf = min^ b* {P,Q + t). In our case, this portion is the construction of As- 

Instead of actually constructing As, we first observe that it suffices to restrict our attention to vertices 
of As, in the sense that each face / of As has a vertex such that all the 1-entries of M(P, Q + f) are 
also 1-entries of M (P, Q -f ^ (the latter matrix can contain additional 1-entries), so it suffices to test for 
reachability in the matrices M{P, Q + C) associated with vertices ^ of As- (Technically, we add to the set 
of vertices one additional point, say the rightmost point, on each disk boundary, to cater to faces that have 
no real vertices.) 

Hence, our parallel implementation of the algorithm will only simulate the construction of the sorted lists 
of vertices along each of the circles Cs{pi — Qj)- Recall that during the parametric searching simulation, we 
collect comparisons that the decision procedure performs and that depend on b, and resolve them. This is 
done by finding the critical values of b at which the outcome of some comparison changes, during a single 
(simulated) parallel step of the algorithm and then by running a binary search through these critical values 
of b, guided by the decision procedure of Theorem |5.6[ In this manner, we maintain a shrinking half-open 
interval I = {a, /3] of values of b that contains bf- Note that we have called the decision procedure at a 
and it has determined that bf > a- Then, as is easily seen, bf must be at least as large as the first critical 
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value of 5 within I (and it cannot be arbitrarily close to a). Assume that we have simulated the construction 
of As, and obtained a half-open interval range I = {a, /?] of <5 that contains < 5 ^. That is, we know that 
a < 6^ < f3, and we know the sorted sequences of vertices of As^ along each circle Cs^{pi — qj). None 
of the comparisons that the decision procedure has performed has a critical value inside I, other than those 
comparisons that have produced (a and) j3. Hence the output representation of As is fixed in the interior 
of I. The rest of the algorithm, which constructs the structure f, traverses the vertex sequences along the 
circles Cs^{pi — Qj), and dynamically updates the reachability data, is purely combinatorial, and does not 
introduce new critical values (i.e., does not involve comparisons that depend on <5^), so there is no need to 
run it at all. Since the decision procedure fails at a and succeeds at /3, it follows that 5^ = /3. 

It is thus sufficient to simulate, at the unknown (5^, an algorithm that 

1. finds fhe infersecfion poinfs of each circle Cs^{pi — qj) wifh fhe circles {Cs:^{pk — Qi) \ Pk ^ P, Qi ^ 
Q}, ofher fhan ilself, and 

2. sorts, for each circle Cs^{pi — qj), the intersection points that were found on its boundary in stepj^ 
along this boundary. 

During the simulation we progressively shrink an interval / = (a, /3] C M that is known to contain <5^. We 
start with I = (0, oo]. 

We first obtain all the 0{m?n‘^) critical values of type[^ sort them, and run an explicit binary search 
among them guided by the decision procedure. (This part requires no parametric simulation.) As a re¬ 
sult I is shrunk to an interval (a, /3], where a, /3 are two consecutive critical values of type|^ This takes 
0{m?n? log(m -|- n) -|- m^n^(l -|- log(n/m)) log(m -|- n)) = 0(m^n^(l -|- log(n/m)) log(m -|- n)) time. 
We can now accomplish step[^ because the property that a pair of circles Cs(p — q), Cs{p' — q') intersect 
either holds for all 5 G (a, /3) or does not hold for any such 5. 

We then execute step|2| The task at hand is to sort, for each circle Cs^{po — qo), the resulting fixed 
sef of infersecfion poinfs along Cs^{po — Qo)- For o^ch pair Cs^{p — q), Cs^{p' — q') of such circles, fhe 
order of fhe infersecfion poinfs can change only af fhe radius 5 of fhe circumcircle of po — qo, p — q, p' — q'. 
We fhen simulate a parallel sorting procedure, fo sorf fhese infersecfion poinfs along Cs^{po — Qo), and run 
if in parallel over all fhese circles. We omif fhe (by now) routine defails of fhis simulafion (see, e.g., [31 
for similar application of paramefric searching). They imply fhaf we can simulafe fhis sorting, for each 
circle Cs^{pi — Qj), using 0{mn) processors and 0(log(mn)) = 0{log{m -f- n)) parallel steps (for a fofal 
of 0{m?'iA) processors). Thus, for each parallel step, we need fo resolve 0{m?v?) comparisons, each of 
which compares (5^ fo a crifical circumradius of fype We run a binary search among fhese crifical values 
using fhe decision procedure. This lakes 0(m^n^(l -|- log(n/m)) log(m -|- n)) time for each parallel step, 
for an overall -|- log(n/m)) log^(m -|- n)) time for 0{log{m?'n?)) = 0(log(m -|- n)) steps. To 

(slightly) improve this running time we use the improvement of Cole Q which finds, for each parallel step, 
the (weighted) median of the (suitably weighted) unresolved critical values involved in this step, and calls 
the decision procedure only at this value, instead of using a complete binary search. This allows us to resolve 
comparisons that contribute at least some fixed fracfion of fhe fofal weigh!, while fhe ofher unresolved crifical 
values are carried over fo fhe nexf step wifh fheir weighfs increased. Proceeding in fhis manner, we make 
only one call fo fhe decision procedure af each parallel step, and add only 0(log(m -|- n)) parallel steps 
fo fhe whole procedure. We fhus obfain an overall algorifhm wifh -|- log(n/m)) log(m -|- n)) 

running lime. 

In conclusion, we gel fhe following main resull of fhe paper. 
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Theorem 6.1. Let P, Q be two sequences of points in of respective sizes m and n, where m < n. Then 
the minimum discrete Frechet distance under translation between P and Q can be computed in 0{m^n‘^{l + 
log(n/m)) log(m + n)) time. 

Discussion. Our algorithm is composed of two main parts. The first part is the construction of the subdivi¬ 
sion As whose complexity is 0{m?n‘^). The challenge here is either to argue that, in favorable situations, 
the actual complexity of As is o{nn?v?), or be able to process only a portion of As that has o{m‘^nf) com¬ 
plexity. Here is a simple illustration of such an approach. Consider the case where P and Q are sampled 
along a pair of c-packed curves, where a curve 7 is c-packed if, for every disk D, the length of 7 n D is at 
most c times the radius of D. Assume also that the sampling is more or less uniform, so that the distance 
between any pair of consecutive points of P or of Q is roughly some fixed value A. We may assume, 
wifhouf loss of generalify fhaf pi = qi. Consider fhe decision procedure wifh a given paramefer 6, and 
observe fhaf if t is any franslafion for which 5*{P, Q + t) <5 fhen ||f|| < 5. Therefore, for each pi G P, 
fhe only poinfs qj G Q fhaf can align wifh pi during a simulfaneous fraversal of P and Q + t, for any 
such “good” franslafion t, are fhose af disfance < 26 from pi. The assumpfions on P and Q imply fhaf fhe 
number of such poinfs is af mosf roughly 26c/ A. Thaf is, insfead of construcfing fhe enfire arrangemenf 
As, it suffices fo consfrucf a coarser arrangemenf, involving only roughly 2c6m/ A disks. Then, fraversing 
fhe coarser arrangemenf is done as before, where each updafe sfep (and fhe following reachabilify query) 
cosf 0(m log(l -|- log(n/m))) fime, assuming fhaf m < n. This improves fhe running fime of fhe decision 
procedure fo O m^(l -|- log(re/m))^, assuming fhaf m < n and 6 < 

Given fhis decision procedure, we can solve fhe opfimizafion problem using paramefric searching. How¬ 
ever, fo ensure fhaf fhe decision procedure does nof become foo expensive, we wanf fo run if only wifh values 
6 = 0(6f). This will become significanf only when 6 f < 60 '■= ^', ofherwise fhe running fime will be 
close fo fhe running fime of fhe algorifhm of Section]^ Therefore, in fhe following we describe how fo solve 
the optimization problem assuming that 6 f < 60 (if the following procedure fails, we run the algorithm 
of Section]^. We also assume, for now, that 6f > A (we explain below how the case where 6f<Ais 
dealt with). We consider the interval (A, 5o) that is assumed to contain 6f, and run an “exponential search” 
through it, calling the decision procedure with the values = 2* • A, for f = 0,1, 2,..., in order, until the 
first time we reach a value 6' = 6i> 6f (and 6' < 6o). Note that the cost of running the decision procedure 
at 6' and at 6f differ by at most a factor of 4, so the cost of running the decision procedure at 6' is asymptoti¬ 
cally the same as at <5^. Moreover, since the running time bounds on the executions of the decision procedure 
at (5i,..., form a geometric sequence, the overall cost of the exponential search is also asymptotically the 
same as the cost of running the decision procedure at <5^. We then run the parametric searching technique 
as above, with the constraint that 6f is at most 6' (i.e., we set 6' as the minimal /3 obtained so far). Hence, 
from now on, each call to the decision procedure made by the parametric searching, will cost no more than 
the cost of calling the decision procedure with 6' (which is asymptotically the same as calling the procedure 

with 6f). We thus obtain an overall algorithm with O m^(l -|- log(n/m)) log(m -|- n)^ running 

time. Note that, in the case where 6f < A, after running the decision procedure with <5 = A, we realize that 
6f < A, and run the parametric searching technique with the constraint that 6f is at most A. In this case, 
the running time of the algorithm is O (c^m^(l -|- log(n/m)) log(m -|- n)). 

The second part of the algorithm presented in this work is the dynamic data structure for maintaining 
reachability in M. It is an open question of independent interest whether this data structure can be improved. 
A related problem is whether the techniques used in our structure can be extended to the general case of 
reachability in planar directed graphs, so as to simplify and improve the efficiency of fhe earlier competing 
mefhod of Diks and Sankowski p8.|. 
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